Mean crossover functions for uniaxial 3D Ising-like systems 
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We give simple expressions for the mean of the max and min bounds of the critical-to-classical 
crossover functions previously calculated [Bagnuls and Bervillier, Phys. Rev. E 65, 066132 (2002)] 
within the massive renormalization scheme of the $ d (n) model in three dimensions (d — 3) and 
scalar order parameter (n = 1) of the Ising-like universality class. Our main motivation is to get 
efficient theoretical expressions to coherently account for many measurements performed in systems 
where the approach to the critical point is limited but yield data which are still reproducible by the 
$d ( n ) model (like in the subclass of one-component fluids). 
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1. INTRODUCTION 

The universal features of the three-dimensional (3D) 
Ising-like systems close to their critical points are now 
well-established by the renormalization group (RG) ap- 
proach jl|. In this theoretical context the Ising-like uni- 
versality is attached to the existence of a unique non- 
trivial fixed point (the Wilson-Fisher fixed point [2] noted 
W-FP in the following) which any Hamiltonian represen- 
tation of an actual system at criticality is driven to under 
the action of the renormalization transformations pj. 

Less known is the existence of theoretical expressions, 
obtained using perturbative field theory (FT) techniques 
\A ti III- lil l&\- which are used to interpolate between 
the critical (non-classical) behavior (controlled by the W- 
FP) and a classical behavior (controlled by the Gaussian 
fixed point, noted G-FP in the following). Such theoreti- 
cal expressions are customarily named classical-to-critical 
crossover functions. 

Actual systems undergoing a second order phase tran- 
sition also display a kind of "classical"-to-critical crossover 
but it is not of the same nature as the theoretical one al- 
luded to above. In actual systems, the "classical" part 
is only a non-critical part (not governed by the G-FP) 
whereas in the theoretical interpolation, the classical part 
still belongs to a critical domain (governed by the G-FP). 
Moreover, the renormalization procedure of FT is reduc- 
tive in the sense that many sources of non-universality 
are discarded because they are unessential in the vicinity 
of the critical point. However, non-universal features be- 
come more and more important as one moves away from 
the critical point. Finally, these non-universal charac- 
teristics are responsible for the distinction between the 
nonasymptotic critical behavior and the asymptotic crit- 
ical behavior. 

Generally, one cannot expect to observe an agreement 



between the theoretical crossover functions and the ex- 
perimental data in a wide domain ranging from the close 
vicinity of the critical point down to a state far away from 
it. Nevertheless, exceptions may exist as in the case of 
the "subclass" of one-component fluids. 

It has been shown [<J, Ojj, uM, LLJ llJ| that, despite their 
wide range of numerical values, the one component fluid 
data collected in the literature regarding a given prop- 
erty (susceptibility, order parameter density, correlation 
length, etc.) can be reduced to a unique scaling curve 
(also called master curve) over a wide critical domain. 
Now, it is well known that this scaling curve can be well 
reproduced by the classical-to-critical crossover functions 
of FT |.l,3f. However, an accurate understanding of the 
range of its validity is lacking. In the present paper, 
we provide modified versions (named "mean crossover 

Suctions") of the theoretical forms given in references 
0, [Tj| . A careful attention is given to the character- 
ization of the Ising like preasymtotic domain where the 
number and the nature of adjustable parameters can be 
easily controlled. A forthcoming paper []j| will illustrate 
the application of the mean crossover functions to provide 
a precise criterion for estimating the effective extension 
of the critical asymptotic domain of the fluid subclass, a 
key point for future developments of a complete equation 
of state (e.o.s.) for fluids. 

Actually, since it is difficult to approach accurately the 
critical point both experimentally and theoretically, the 
effective (observed) scaling domain is limited towards the 
critical point whereas, in the opposite direction, the scal- 
ing curve stops when the systems are no longer critical 
and before they display any kind of "classical behavior". 
Here classical behavior means the "critical" behavior con- 
trolled by the G-FP. Consequently, most of the useful 
data on a system near criticality belong to an intermedi- 
ate (nonasymptotic) regime where it is not valid to repro- 
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duce the data with the currently used expansion limited 
to a pure power law eventually corrected by one Weg- 
ner's term [lfl] (the validity of such a form implies the 
close vicinity of the critical point). 

For development of an e.o.s., it is of great importance 
to have theoretical expressions which reproduce the phe- 
nomenologically observed scaling curves and which allow 
one to situate unambiguously any data set with respect 
to the actual (system dependent) field distance to the 
critical point. It must be stressed that, in this view, the 
aim is not to test the most refined theoretical hypotheses, 
or the most precise estimates of the universal numbers as- 
sociated with the asymptotic or pre-asymptotic critical 
behavior. In addition, our scheme only applies to the 
primary critical paths along the critical isochore in the 
homogeneous and the non homogeneous domain, with- 
out any supplementary consideration on the roles played 
by other nonuniversal scales associated with the eventual 
effects such as non symmetry pj| and mixing [11] of the 
scaling field variables very close to the critical point. 

>From a pioneering study []jj] of the available data on 
xenon [2fl W& > ^ is known that the classical-to-critical 
crossover functions obtained from the <I>^ (n) model (with 
d = 3 and n — 1) fit the one-component fluid results, in- 
troducing explicit determination of the fluid-dependent 
scale factor for the relative temperature field. More re- 
cently, a similar study J2£ of several measurements near 
the liquid-gas critical point of 3 He J23, has provided an 
equivalent conclusion, demonstrating that, 3 He as a sim- 
ple fluid, is also situated very close to a renormalized 
trajectory that links the G-FP to the W-FP. However, 
the authors of Ref (j. y. [fj. |lj| have given a practical 
(max or min) form of their effective functions which is 
better adapted to the eventual experimental test of the 
"best" theoretical estimates of the asymptotic universal 
quantities (attached to the close vicinity of the critical 
point), rather than to an unambiguous characterization 
of a nonasymptotic critical domain of a given subclass of 
systems. In fact, accounting for the error estimates at- 
tached to the universal critical quantities under two max 
and min bounded sets of functions prevents the charac- 
terization of the scaling curve clearly which, of course, 
is currently accessible experimentally at some non-small 
distance (to be determined) to the critical point but also 
is sometimes accompanied by a relatively poor accuracy 
in the measurement. Moreover, due to the decrease of 
the theoretical error as the distance to the critical point 
increases, the mean values of the bounded functions are 
sufficient to characterize the scaling curves. 

Our aim in the present paper is to determine in- 
between controlled functions that are better adapted 
for unambiguous determination (within a current experi- 
mental mean accuracy) of the characteristics of the scal- 
ing curves of the subclass of one-component fluids. An 
important point to note a posteriori is the unchanged 
value of the temperature scale factor determined in [lflj . 
However the uncertainty is considerably reduced showing 
that it was essentially generated by the theoretical uncer- 



tainties of the universal values calculated in the asymp- 
totic regime. That demonstrates the need to provide 
mean classical-to-critical crossover functions which, by 
construction, have a well defined single asymptotic limit 
(i.e. "Ising-like-well- defined") . We will then be able to 
determine significative values of the scale factors which 
characterizes the scaling behavior occuring in the inter- 
mediate critical domain. 

The theoretical crossover functions which we are inter- 
ested in have been derived from a massive renormaliza- 
tion (MR) scheme applied to the <b A d (n) model. The MR 
scheme has been initially developed in references 0. LJ Lj| . 
hereafter referenced MR6, using the Borel resummation 
technique based on the results of the sixth-loop series 
[2J|- They have been recently revisited in reference [14| . 
hereafter referenced MR7, to account for an extension 
to the seventh-loop series [2||. We do not consider here 
the crossover functions determined by Dohm and cowork- 
ers 0, ll| who have used another renormalization scheme 
(minimal) within which the known series are shorter than 
in the MR scheme and have provided their crossover func- 
tions under implicit forms. 

The paper is organized as follows. 

In Section 2, we introduce the main characteristics of 
the MR scheme and of the crossover functions. Special 
attention is given to the non-universal nature of the ad- 
justable parameters that are introduced in the max and 
min theoretical functions for fitting experimental data. 

In Section 3, the mean functions are determined rely- 
ing on the properties of the theoretical functions in the 
two limiting 3D Ising-like and mean field-like descrip- 
tions, respectively close to the W-FP and G-FP. Such 
descriptions correspond to the pre-asymptotic domains 
(PAD) near each fixed point where a Wegner expansion 
restricted to two terms (leading and first confluent terms) 
is valid. The Ising-like PAD includes the correlations 
between parameters due to the error-bar determination 
of the exponents and amplitude combinations very close 
to the W-FP. The addition of the equivalent mean field 
PAD description very close to the G-FP leads to a mean 
crossover function with a limited number of calculated 
parameters for this function (three in the selected ex- 
ample). A well-controlled form of any mean crossover 
function can be obtained in a similar manner, i.e. three 
calculated parameters for each function. In such a sit- 
uation, the theoretical crossover forms are obtained for 
a unique value of one parameter among the three. This 
parameter acts then as a relative sensor to estimate the 
dominant nature, either (Ising-like) critical, or (mean 
field-like) classical, of the calculated crossover. Using 
this sensor, we propose an explicit criterion to measure 
the extension of the Ising-like PAD along the critical iso- 
chore where a four parameter characterization of each 
system is well-understood in scaling nature. In section 4, 
we provide a conclusion. 

In Appendix A, assuming knowledge of the criti- 
cal temperature of the system, we illustrate the three- 
adjustable-parameter characterization of the comparison 



with experimental data. The emphasis is on the role of a 
characteristic microscopic length scale, the reminiscence 
of which is carried by the dimension of the "bare" (f) A - 
coupling go [see Eq. (QJ below]. 

Appendix B gives some details on the derivation of the 
mean crossover function for the particular case of the or- 
der parameter in the heterogeneous domain, with a view 
to better account for large (theoretical and experimental) 
uncertainties in the determination of its corresponding 
first confluent amplitude. 



2. THE MASSIVE RENORMALIZATION 
SCHEME 

2.1. The model 

The calculations of the crossover functions in the <&^ 
model rely upon the renormalization program of per- 
turbation FT. This scheme makes the perturbative ex- 
pansion of the correlation functions free of ultra-violet 
divergencies. Before renormalization, the dimensionless 
Hamiltonian (a true hamiltonian divided by fc^T, with 
ks the Boltzman constant) of the actual system reads: 
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in which d is the space dimension and h and $o are vec- 
tors of dimension n. h is the (magnetic-like) ordering 
field vector. $o is the spin-like vector. The ordering 
field vector h is not renormalized and so will no longer 
be considered explicitly in the following (except when 
required). The coupling field r$ and the coupling con- 
stant go are the bare (physical) parameters, which are 
two system-dependent quantities characterizing the crit- 
ical point location of the physical system. 

For the sake of simplicity, we limit ourselves in the 
following discussion to the scalar case n = 1. 

At this stage, it is worthwhile to state that the dimen- 
sion of any quantity appearing in TC is expressed in terms 
of only one inverse length unit: a wave- vector cutoff, A 
(or the inverse of a lattice spacing for magnetic solid sys- 
tems for example). Since H is dimensionless, a simple 
evaluation of Eq. JQ shows that the dimension (in unit 
A) of $o is =j^, while that of r is 2, and that of go is 
4-d. 

When d = 4, the perturbative expansion in powers 
of the dimensionless coupling constant go involves ultra- 
violet divergencies which are removed by redefining the 
three initial Hamiltonian (bare) parameters {$o,^o,ffo} 
into new (renormalized) ones {$,m, g}, via the following 
relations: 



$ = [Z 3 (g)]*$ 



,— -,4-d Z\ (g) 
9a = (m) g- 
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with subtraction conditions which, in four dimensions 
(d = 4), are required to make the renormalized pertur- 
bative expansion in powers of g finite. 

In three dimensions, however, only the mass renormal- 
ization is needed to have a finite theory. Consequently, 
with the introduction of rn instead of ro, the renormaliza- 
tion functions, like Z\ and Z%, even expressed in terms of 
the "bare" coupling go and without ultra-violet regulator 
A (set to infinity), are well defined functions. Introduc- 
ing the notation T Q ' ({q,p};m, g ,d) for the Fourier 
transforms of the bare Appoint vertex functions with L 
"insertions" of the squared field <&q, the following defini- 
tions stand (for d < 4) 
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Z^=T^>({0};m,go,d) 

in which Zi is a renormalization function which restores 
the linear measure of the distance to the critical temper- 
ature T c , originally defined by the bare parameter ro and 
which has been lost when introducing the renormalized 
mass m [via Eq. (|2j and a condition, not written here, 
which defines rn as the inverse correlation length]. 

Since A has been eliminated in the renormalization 
process, m plays the role of the effective wave- vector scale 
of reference and could be used as the unit to express the 
dimension of any dimensioned quantity. However, be- 
cause rn (the inverse of the correlation length) vanishes 
at the critical point it is preferable to use go which, for 
d 7^ 4, is a dimensioned constant at the critical tem- 
perature [see Eq. (1 1 01) below]. Hence, (go) 3 ^ 3 will be 
substituted to A to play the role of the wave- vector unit 
(see below). Notice that in three dimensions go has ex- 
actly the dimension of A (i.e. of the inverse of a length) 
which is very convenient but not essential. 

The dimensionless renormalization functions Zi have 
been calculated u p to sixth order [2J| and then partly up 
to seventh order J2a| in powers of g. These series have 
been summed to estimate the critical exponents with 
great accuracy [2a, |23 and also to determine nonasymp- 
totic critical functions in the homogeneous phase_p,0|. 
With the calculation of supplementary integrals jfj, the 
calculations have been exended to the inhomogeneous 
phase up to fifth order allowing the determination of the 
nonasymptotic critical functions 0] in this phase and of 
the equation of state [271 12^ |. The calculations of the 



nonasymptotic critical functions have been revisited [Tj] 
in order to account for the most recent estimates of the 
asymptotic universal critical quantities [23 and also to 
provide complete classical-to-critical crossover forms [l3 | 
which we are presently interested in. 



2.2. The crossover functions 

In the perturbative framework, the renormalized cou- 
pling g may take on any value in the range [0, g*] where 
g* is its W-FP value (a value which may be estimated 
in three dimensions by looking at the nontrivial zero of 
some series) and is its G-FP value. The crossover func- 
tions have been obtained by resumming the series of a 
physical property P of interest such as the correlation 
length £(g) [— (m)~ ] or the susceptibily x(g), for a 
discretized variation of g in this range. An expression 
for the critical scaling field t (g) (which "measures" the 
physical distance to the critical point) is obtained via an 
integration of the resummed series of Z 2 (g) (see ref. Q 
for example). Hence the variations of £ or x, in terms of 
t, are primarily obtained implicitly via the dummy pa- 
rameter g. Explicit functions of t representing £ or \ are 
then obtained by fitting ad hoc forms to their discretized 
evolutions in the range ]0, g*[. 

It is to be noted that, appart from r$ which is used to 
determine t (see below), the only remaining dimensioned 
parameter of the $^ model is the bare coupling go with a 
dimension (in length unit) equal to d — 4 . Then, all the 
final functions are reduced by the appropriate powers of 
<7o so as to be dimensionless. 

In MR6, contrary to MR7, the entire crossover (cor- 
responding to the complete range ]0,g*[) had not been 
published. Also, in MR6, the error analysis was made 
independently from that associated with the estimates 
of the critical exponents done by Le Guillou and Zinn- 
Justin 26]. In MR7 the convergence criteria for the Borel 
resummation of the different functions, like £ or \, was 
chosen such that the max and min bounds of the re- 
sulting critical exponents agreed as closely as possible to 
the (revised) values obtained by Guida and Zinn-Justin 
[23- However, in doing so, it is more than likely that the 
error bars have been over estimated compared to what 
the resummation method used would have naturally indi- 
cated (following the rules applied in MR6). This is why 
the results of MR6 are also of interest to us especially 
in the case of the order parameter in the heterogeneous 
domain for which the large uncertainties of MR7 are not 
favourable for accurate fitting of the experimental data 
(see appendix B). 

In MR7, two tables were presented to give an envelope 
for each function Fp (if) representing a (dimensionless) 
model property P versus a discretized dimensionless scal- 
ing field t* [defined in Eq. l(TT|l below] over the entire 
range ]0, g*[ for the max and min bounds. The follow- 
ing expressions Fp (t*) were used to continuously fit the 



discrete data from each particular table 
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with 



D$ (t*) = A - 1 



S± VF+1 
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In Eq. (2J, 3 < K < 5, depending on the required fit 

quality for each property P = I (x*)~ ; (£*)~ ;C*;m*> 

[the superscript * indicates dimensionless quantities, x* 
is the susceptibility; £* is the correlation length; C* is 
the heat capacity; m* is the order parameter in the non- 
homogeneous domain (note the distinction with the dec- 
orated to used for the renormalized mass)], ep and A 
are the leading and first confluent universal exponents, 
respectively. The symbol ± indicates the possible homo- 
geneous (+) and nonhomogeneous (— ) domains. All the 



constants Z D , X 



,,. > ,. j, Y Pi , and S Pi in Eqs. © and © are 



tabulated in [Tj]- They result from the fitting of eqs. 10 
EJ to the theoretical calculations done point by point in 
the discretized complete t* = {oo,0} range correspond- 
ing to the complete range ]0,g*[. The critical behavior 
of the specific heat is particular such that, compared to 
Eq. (§J, it involves an additive critical constant that we 
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2.3. Physical validity of the functions 

2.3.1. Analytical corrections discarded 

To compare the theoretical functions to measurements, 
we must relate the Hamiltonian parameters (ro, go, h) to 
their physical counterparts. This is done by the (usual) 
basic assumption that the bare quantities (go, ro, h) are 
analytical functions of the corresponding physical quan- 
tities T and H (in usual notations for magnetic systems) 
which control the approach to the actual critical point. 
Implicitly, we admit that the Hamiltonian energy is com- 
parable to the physical free energy measured in unit of 
kpT = kpT c very close to the actual critical temperature 
T c of the system. 

Assuming that analytical corrections to scaling are 

negligible, this introduces two arbitrary scale factors 

(noted i? and tp in the following) associated with the two 

bare fields ro and h, respectively, and one constant in- 

i 
verse length scale (go) i ~ d which fixes the dimensionality 

of the physical variables. This can be accomplished as 

follows: 

1. At h = 0, the bare field ro = ro (T) must be related 
to the actual temperature T of the system, leading 
to define its critical value r 0c = r (T c ) from the 
actual critical temperature T c of the system. For T 



close to T c , it becomes possible to relate the bare 
field difference r (T) — r 0c (T c ) to the temperature 
distance T — T c by the following analytical (linear) 
approximation 



ra - r Qc = const x (T - T c ) + O (T - T c ) 



(7) 



In an equivalent manner, at rg 



roc 



the bare field 



h must be related to the actual (magnetic field) 
variable H (in notations for magnetic system) of 
the system, leading also to the possible analytical 
(linear) approximation 



h = const xH + 0[H 2 ] 



(8) 



for h close to zero. Correlatively, the linearization 
between the bare order parameter to = ($o) to 
the actual (magnetization) density variable M (in 
notations for magnetic system) reads 



m= constxM [(T-T C ),H] + 

O \{T - T c f ,H 2 ,H x (T - T c 



(9) 



where the two constant prefactors of Eqs. © and 
© are interrelated by thermodynamic considera- 
tions attached to the conjugated variables M and 
H (or to and h equivalently) . 

Finally, at the above linearized order of the rela- 
tions between the bare fields and the physical fields, 
close to the critical point defined by T — T c = (or 
r o — r 0c = 0) and H = (or h = 0), we can intro- 
duce one finite value of the coupling constant such 
as, 



go = const + O (T - T c ) + O [H] 



(10) 



which appears then as a system-dependent critical 
quantity which must take the d — 4 dimension (in 
length unit). In the Eqs. (0 to ©, one generally 
keeps only the leading terms but one must keep in 
mind that far away from T c , the neglected second 

order analytical terms O (T — T c ) , O [i? 2 ] , and 

O [H x (T — T c )] (within ro, h and to) could have 
some importance especially in approaching the G- 
FP. 

With go defined by Eq. 11 Oil , we can now appropri- 
ately make dimensionless the bare quantity ro by 
introducing the dimensionless scaling field t* used 
in Eqs. (0 and ©: 



t* = 



ro - r Qc 



(.90 



(11) 



For the actual system, it is convenient to define the 
reduced temperature distance in units of the criti- 
cal temperature (already chosen to express the en- 
ergy unit of the hamiltonian) , leading to the usual 



._ T-T c 
T 

and At 



As- 
are 



notation of the thermal like field At* 

suming a "sufficiently" small At* , t * 

related by the first dimensionless arbitrary scale 

factor ■&: 



t* = -d x At* 



(12) 



Similar assumptions stand for the second scaling 
field h which must be related to H by Eq. ©• 

A second arbitrary scale factor appears, noted tp, 

hrr- and 



between the dimensionless field h* = - u 

its corresponding dimensionless physical quantity 
H* (here in notations for magnetic systems): 



4;xH* 



(13) 



Correspondingly, the dimensionless order parame- 
ter field m* = ($*,) = 



™_ 2 is then related 

(go) 2T4=dT 



to the dimensionless physical quantity M* (here in 
notation for magnetic systems) by 



to* = (ipy 1 x m* 



(14) 



Finally the "length" (.go) 3 ^ must be related to a 
microscopic length scale a characteristic of the ac- 
tual system as: 



(go) 11 ^ =u*xa 



(15) 



where Uq is a dimensionless number similar to ■& 
and tp. Notice that Up has the characteristic of tak- 
ing a value which depends on the (a priori) choice 
of a system dependent length scale a, not attached 
to the critical behavior. That characteristic will 
be used to determine a subclass of "comparable" 
systems [i.e. systems having a comparable charac- 
teristic length scale a (see below)]. 

As will be shown in Appendix A, it is worthwhile 
already indicating here that the consideration of a 
via Eq. (1 1 5 1 is not required when the correlation 
length is considered alone. This is because the di- 
mensionless correlation length £* of the theoretical 
model is naturally compared with the experimental 
measurement £ via the following relation: 



e* = (9o) T = 3 *Z 



(16) 



In Appendix A, limiting our calculations to the d = 3 
case, we illustrate the way the adjustable parameters are 
determined where the emphasis is also on the choice of 
the microscopic length scale a appearing in Eq. ltl5ll . 
When the critical temperature of the system is known, 
fitting the asymptotic two-term expansion of our theo- 
retical crossover functions to the singular behavior of the 
experimental quantities, permits unambiguous determi- 
nations of the three system-dependent parameters $, tp, 
and go- We then indicate how the dimensionless coupling 



constant Mq may be used to characterize "comparable" 
systems which belong to a same subclass of universality 
when an explicit length scale unit (a) is known for each 
system. We show that the theoretical crossover functions 
then provide an explicit analytical form of the "master" 
(i.e. unique) singular behavior associated to this subclass 
and discriminate the role of the energy and length scale 
factors for the subclass. 



2.3.2. Non-analytical corrections discarded - 
Some ideas on renormalization 

The renormalized $^ (n) model of the perturbative FT 
is an efficient reduction of a complicated mathemati- 
cal problem which originally involves an infinite num- 
ber of parameters to one parameter g (the renormalized 
<& 4 coupling g of section l2~T|) . In order to better under- 
stand the impact of this reduction, it is necessary to con- 
sider the complete and non-perturbative renormalization 
group theory developed by Wilson Q. 

Explaining the use of the renormalization theory in 
the study of critical phenomena is out of the scope of 
the present paper. However, it is worthwhile to indicate 
briefly some general ideas which may help someone to un- 
derstand the use of our crossover functions in comparison 
with experimental data. 

The RG theory is a general theory to treat situations 
where infinitely many degrees of freedom are correlated, 
such as near a critical point where the correlation length 
£ diverges. Since in such cases, we are essentially in- 
terested in describing the large distance behavior, and 
because of the correlations, one may represent the state 
of a system near criticality by means of a dimensionless 
hamiltonian TC which depends only on a local field <fi (x) 
which summarizes, over a volume of linear size A -1 cen- 
tered on x, only general properties of the genuine local 
variable which critically fluctuates (e.g. the spin vari- 
ables Si of the Ising system). Beyond 4>{x), TC depends 
also on the relevant physical parameters like the temper- 
ature T and the magnetic field H (for the sake of sim- 
plicity we shall not consider H) . The form of TC is quite 
general provided it satisfies some required properties of 
symmetry. For systems which are 0(l)-symmetric [more 
precisely ^-symmetric, i.e. invariant under the change 
4> — > —(f) such as Ising- like systems], then TC must be an 
even function of <j> (when H = 0) . 

In order to concretize a bit the form of TC, one may 
look at its expansion for small values of <f> {x) , in the case 
of the ^-symmetry, it becomes (do being a term which 
is usually discarded in FT) : 

TC = a a + d d x |oi (V0) 2 + a 2 4> 2 + a 3 4 

+a 4 </» 6 +a50 2 (V0) 2 + ---} (17) 

in which the expansion in powers of the derivatives is 
a consequence of the short range interactions between 



the original spins and the coefficients a, depend on T. 
Actually one is free to redefine the global normalization 
of the field so as to set oi = constant (in general one 
chooses a\ — |) but this is not mandatory. Furthermore 
all the dimensions are measured in terms of A, so that it 
is convenient to deal with dimensionless quantities in TC. 

The set < ai ( if- J > characterizes a given physical system 

near a critical temperature T c . It may be seen as the 
coordinate of a point in a space S of infinite dimension. 
The problem is that calculation of any critical quantity 
with TC is very complicated. The reason is that the pa- 
rameters at hand < ai ( if- J > are attached to the length 

scale A -1 , while the physics under consideration refers 
to phenomena that occur at all the length scales smaller 
than £ and greater than A -1 . Now because the ratio 
^4t tends to infinity, there are an infinity of scales to be 
accounted for to solve the problem. 

Introducing a parameter s such as s £ [0, oof, the 
renormalization group transformations precisely con- 
struct successive hamiltonians Ti(s), obtained by inte- 
grating out the degrees of freedom over the scales ranging 

in (e~ s A) , A -1 , and then rescaling back the wave- 
vector scale (hence A' = e~ s A — » A). By construction, 
the renormalized hamiltonian TC (s) presents a correla- 
tion length £ s which is reduced, compared to the £ of 
the initial hamiltonian TC (s = 0) = TC, with £ s = e~ s £, 
and the effective hamiltonian TC (s) is no longer critical 
when j=t ~ 1. If the initial hamiltonian is critical (i.e. 
if £ = oo), then TC (s) reaches a fixed point when s — > oo. 

The evolution of TC (s) under an infinitesimal change 
of s is governed by an equation, called the exact RG 
equation. It is a complicated integro-differential equation 
which may be expanded in powers of the derivative of the 
field without losing the non-perturbative character of the 
complete theory (for a review see [23]). 

At leading order of the derivative expansion, called the 
local potential approximation (LPA) , the exact RG equa- 
tion reduces to an ordinary differential equation for a 
function V of one variable <fi. Its relation to the $^ model 
is better seen if one expands V in powers of <j>. Imposing 
the property of parity in <f>, it becomes (for <f) close to 0) : 

V (0, s) = V (a) + V 2 (s) 2 + Vi (s) 4 + V 6 (s) cj) 6 + ■ ■ ■ 

(18) 
The form of Eq. I|18|l shows the relation with Eq. QJ 
when the derivatives of the field are neglected (again, for 
the sake of simplicity, the account of h, which would re- 
quire the presence of odd powers of <f), is not considered 
here). In particular, the initial form chosen for V may be 
directly related to the bare parameters of the $\ model 
with, e.g. A 2 V2(0) = f 3 and A 4 -^ (0) = § (and 
Vo (0) = ^6 (0) = 0), they are parameters attached to the 
microscopic length scale A" 1 of the initial hamiltonian 
(s = 0). As for the renormalized coupling g of the per- 
turbative framework, it is not a constant but a running 
parameter like V4 (s). To be more precise, one must look 
at the flow of V {<fi, s) solution of the RG equation. 



In LPA, the derivatives of the field are neglected but 
one may account for all the effects of all the Hamiltonian 
terms via the evolution of the complete function V ((/)) 
under the action of the RG transformation. 

In three dimensions, a non-trivial fixed point V (</>) ex- 
ists and may be (numerically) determined. One may il- 
lustrate this by writing (as a formal polynomial of (j>): 



V {<t>) = V + V2<t> 2 + Vi<t> A + V e cf> 6 + 



(19) 



and then calculating the values of Vq, V2, V4, Vq, etc. 
As a consequence, one may situate the fixed point in 
the formal space S truncated to the set of couplings 
{Vq, V2, V4, Vq, ■ ■ ■ } of infinite dimension. To visualize the 
evolution of a RG trajectory of V (</>, s) (the evolution of 
V (<p) under a continuous RG transformation controled 
by the parameter s <E [0, 00 [ in S), one may project it in 
the plane {V4,Ve}, for example. 

We may choose to integrate the differential equation 
with an initial simple form for V (4>, s — 0) like, for 
example: 



V (0, 0) = V 2 (0) <\? + V 4 (0) ^ + V 6 (0) </> 6 



(20) 



in which V2 (0), V4 (0) and Ve (0) are numbers represent- 
ing the initial values of V (s). As soon as s ^ 0, then 
V {<p, s) contains all powers of </>. 

To approach the fixed point V (4>) starting with (1201 , 
one must adjust one of the three initial parameters, e.g. 
p2(0), to a particular value V<2 C , which depends on the 
value chosen for the other parameters 14(0) and Vq(Q). 
The requirement is similar to the adjustment of the crit- 
ical temperature with a view to reach the critical point 
in actual systems. The resulting value V c ((f)) is said to 
belong to the critical subspace S c of S, which forms the 
domain of attraction to the W-FP. S c has the same di- 
mension as S minus one (the codimension of S c is equal 
to one) which corresponds to the fact that the W-FP 
has only one direction of unstability, corresponding to a 
direction locally orthogonal to S c (we suppose h = 0). 

It is thus possible to illustrate the RG trajectories in 
S c which, starting from some arbitrary points of S reach 
the W-FP when s — > 00 (one initial parameter must be 
finely adjusted). This is demonstrated in Figure 1. On 
this figure one sees that, whatever their starting points, 
the critical trajectories reach the W-FP asymptotically 
along a unique ideal trajectory which links the G-FP to 
the W-FP. The closer the starting points are chosen to 
the G-FP, the longer is the way along this ideal trajec- 
tory. This ideal RG trajectory is called the "renormal- 
ized trajectory" (RT) because the RG flow running on 
it corresponds, in the vicinity of the G-FP, to the RG 
flow discovered in the perturbative theory of renormal- 
ization. Since the RT is a manifold of dimension one, a 
single parameter (the renormalized coupling "constant" 
g) is sufficient to characterize the RG flow running on it. 

Hence, in some sense (there is an arbitrariness in the 
definition of the renormalized coupling) the renormalized 
(or running) coupling g of the perturbative framework 




Figure 1: Ideal Renormalization Group trajectory in the bi- 
nary diagram of coupling parameters V4 and Ve [see text and 
Eqs. (EJ and 12TB ]. 



implicitly "flows" along the RT. Actually, g is a function 
of s, more precisely in the circumstances of the massive 
field theory, it is a function of ^ . The difference with the 
complete Wilson's RG theory is that, in the perturbative 
framework, one has implicitly assumed to be "exactly on" 
the RT, so that the distance to the W-FP is measured by 
the values taken by a unique parameter g (^) — g* (g* 
being the value of g at the W-FP). The consequence is the 
presence of only one family of correction-to-scaling terms, 
such as in the following expansion of the (dimensioless) 
correlation length for example: 



t(t*) =^ ± |**| _ 



E 6 " 
fc=i 



t*\ 



/.■A 



(21) 



with A a specific exponent the value of which (close to 
0.5) may be estimated by resumming series of the renor- 
malized perturbation theory of the $^ model. It is pre- 
cisely expansions like in Eq. ll2l|l that have been summed 
into the crossover functions considered above. 

The other transients not accounted for by the $^ model 
(by implicitly assuming that the only possible RG tra- 
jectory is the RT) are responsible for other kinds of 
correction-to-scaling terms characterized by a hierachi- 
cal set of exponents Aj (j = 2, • • • , 00) such that: 

A < A 2 < A 3 < ■ • ■ < Aoo 



The fact that A (the exponent which controls the 
asymptotic approach, along the RT, to the W-FP) is the 
smallest exponent is well illustrated on Figure Q] by the 
coincidence of all trajectories with the RT before reaching 
the W-FR 

The few estimates of A2 indicate that its value in three 
dimensions and for n = 1 is of order 2 A so that, on a 
general ground, the expansion IJ21H is a priori not valid 
beyond the first correction term. This is why one often 
uses the following two-term Wegner expansion to analyse 
the experimental data: 



£(At*) 



£ l A ^* 



1 



- af I At* 



(22) 



with the exponents v and A fixed to their FT values. 
The domain of validity of such a two-term expansion is 
called the preasymptotic domain (PAD) and noted |i*| < 

plsing 
'-'PAD- 

For the same reason, the crossover functions calcu- 
lated in FT were not expected to generally reproduce 
the experimental data beyond the Ising like PAD. How- 
ever this objection does not account for the amplitudes 
of the correction terms. It may well occur that the co- 
efficients of the corrections associated to the exponents 
Aj (j = 2, • • • , 00) are very small so that in effect only 
the <I>^-like corrections have to be considered even on a 
range of t* for which several (more than two) terms of the 
expansion of Eq. 1I2D1 are not negligible. It is the case of 
initial Hamiltonians the coordinates of the critical point 
of which correspond to a point lying very close to the RT. 
This particularity may even validate the entire crossover 
functions of FT if the initial point lies close to the G-FP. 

Of course, the RT does not represent the unique pos- 
sibility of approaching to the W-FP and the infinite di- 
mension of S c leaves room for an infinity of possibilities. 
For example, there exists another attractive trajectory 
which is characterized by an asymptotic approach to the 
W-FP controlled by A2, and supplementary adjustment 
of the initial Hamiltonian is needed to obtain this kind 
of approach. For systems represented by such Hamilto- 
nians, the crossover functions of FT have no utility at 
all since none of the correction-to-scaling terms is cor- 
rect. Similar and more common are those systems which 
correspond to an approach to the W-FP from the "op- 
posite side" compared to the RT. Their critical behavior 
are characterized by amplitudes of the first correction- 
to-scaling (controlled by A) with a sign opposite to that 
generated by the approach along the RT 30]. 

Up to now, it is not possible to say a priori which kind 
of approach to the W-FP may correspond to an actual 
system which belongs to a given universality class. A way 
to find out is to try a fit of test functions to experimental 
data. From such tests, it seems that the subclass of one- 
component fluids corresponds to Hamiltonian lying very 
close to the RT [31J- This is why the crossover functions 
of FT are good candidates to fit the scaling curves of that 
subclass. 



3. THE THEORETICAL MEAN CROSSOVER 
FUNCTIONS 

In this section, we determine mean values of the 
crossover functions. We also take the opportunity of sim- 
plifying the ad hoc functions by limiting the products in 
Eq. (0 to only three terms. Hence we have to determine 
the value of parameters like those entering Eqs. (© and 
© , but with K = 3 for all quantities and so that the cor- 
responding ad hoc functions lie just in-between the max 
and min bounded functions published in Ref. [1J|- To 
accomplish this, the following is taken into consideration. 

1. The required precision in fitting each property re- 
quires the recourse to phenomenological confluent 
functions D^ (t*) [Eq.©] which essentially account 
for the crossover between the values A ~ | of the 

0, and the 

1 of the mean-field-like "conflu- 

-> 00. Such a crossover intro- 



Ising-like confluent exponent as t 

exact value A. m f = \ 

ent" exponent as t* - 

duces the following condition 



S$,i = S$,2 (1 - A + A m/ ) 



(23) 



which confers to one parameter among S P 1 and 
S P2 , a noticeable difference from the other con- 
stants appearing in Eqs. JU and ©■ 

2. In [14| . the ad hoc functions for the susceptibility in 
the homogeneous phase are already defined with a 
three-term product. In that case, the behaviors of 
D+ (£*) as a function of t* as illustrated in Figure 
clearly show that the condition 



« = 1 






can be used as an indicative sensor of the classical- 
to-critical crossover (C3) domain (£*, stands for the 
values of t* where the ad hoc confluent functions 

take on precisely the mean value Ai = j 22 ^)- 

The extension SCcz of the associated "intermedi- 
ate" t* -range can be measured by a characteristic 
amplitude w < 1 [which remains to be defined, 
see below Eq. IJ391) ]. such that to m (S^ 2 ) < 



t* e [5C C3 ] < i (S* 



£ (5± 2 ) 2 (see Figure©. 



j 2 j , separating thus unam- 
biguously an Ising-like asymptotic domain of ex- 
tension t* < £ Ism 9 = w (S^) and a mean-field 
like asymptotic domain of extension t* > C mf — 

'x.: 



3. The leading amplitudes Z p , associated with their 
respective universal exponents ep satisfying scal- 
ing laws, are unambiguously related by universal 
amplitude combinations between them. 

4. A part of {Xt, l^j-p's accounts for the universal fea- 
tures associated to the critical confluent corrections 
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Figure 2: Typical crossover behaviors of the confluent 
crossover function Dj (£*) [Eq. J^J], for the bounded (max 
and min) critical confluent exponents. as a function of 
the thermal field (at zero ordering field) for the suscepti- 
bility in the homogenous phase (T > T c ). The vertical 
lines are indicative of typical discrete variations of the con- 
fluent function, showing the practical interest of the condi- 
tion £+2 (A, A m /) y^ = 1 (crosses), which corresponds to 

Z)J (to) = Ai = ^r 1 ^-- Using a single perfactor vj < 1 

which remains to be estimated [see below Eq. I|4()fl]. it is 
possible to define the extensions t* < £, Iatna — zu (S^ 2 ) i 
t* > £ mf = i (StoV 2 , and C Isins < t* G [SC C3 ] < C mf , of 



X-2J 



either, the Ising-like (I sing), or the mean-field- like (mf), or 
the classical-to-critical crossover (C3) domains, respectively 
(see text for details). 



to scaling, explicitly characterized by the single uni- 
versal exponent A and one first confluent amplitude 
(see below) for t* — ► 0. 

5. Another part accounts for the asymptotic mean- 
field behavior for t* — > oo [implicitly characterized 
by the two (leading and confluent) exponents ep, m f 
and A m f and their associated amplitudes (see be- 
low)]. 

6. The remaining part accounts for the expected 
classical-to-critical crossover (C3) in the interme- 
diate i*-range (which also remains to be defined, 
see point 2 just above). 

These remarks provide constraints on the parameter 
entering each specific ad hoc functions. Let us consider 
those constraints explicitly. 



3.1. 3D Ising-like PAD description 



Ising-like PAD descriptions given by the bounded func- 
tions F P<max and Fp iTnin , respectively. 

Starting with Eq. ©, we first consider the two-term 
(Wegner) expansion of Fp (t*): 



pPAD /,* 
r PJsing \ L 



z$(t* 



1 



Z^ ± (i*) A 



(24) 



where the function Fpf? n (£*) is valid within some 

Ising-like PAD extension t* < £ p s ^. To estimate mean 
exponents and mean amplitudes of Eq. IJ24H from max 
and min ones, we impose the following four obvious con- 
ditions 



e P 



ep. 



e P , 



7L P — yZ PmM Zp min 



A = 



A, 



7l,± 



V^3 y-zt -y± 

~ l^i=l ^ L P,i 1 P,l 



,+xl 



(25) 



(26) 



(27) 



(28) 



The leading term of our mean function is then identical 
to the leading term of the mixing function 



F E 

P. mix 



{Fp, 

{Fpn 



t*Di 



1-E 



^'Dpmtx (**) - j Fp t max t*,Dp max (**) j 

, (<*)]} 

(29) 
proposed in [lj], with E = |. The three first Eqs. 
II25H . 12(ili . and IJ27H . provide unequivocal determination 
of the mean values of the three parameters ep, Z p , and 
A, respectively. The Eq. l|28jl involves exclusively the 



3.2. Mean-field-like PAD description 

To recover the (asymptotic) mean-field-like behavior 
of Fp [t*,Dp (**)] in the limit t* -» oo, the following 
relations are required 



ep, m f = ep 



2 2^, Y P,i 



(30) 



Since the error estimates are the largest in the vicinity 
of the W-FP, it is in the Ising-like PAD limit that the 
determination of the mean value has important conse- 
quences. 

Let us define a set of four constraints to control the 
mean Ising-like PAD description within the two limiting 






(31) 



where ep t mf and Zp . are the mean-field (classical) val- 
ues of the asymptotic exponent and amplitude, respec- 
tively. Such a mean-field situation as t* — > oo, can be 



10 



easily characterized in an equivalent manner to the above 
Ising situation as t* — > 0. 

Therefore, the following restricted two-term expan- 
sions 



T?PAD 



(t* 



z 



P,mf 



(t* 



l + Z 



1,± 

P,mf 



(t*y 



(32) 
can be easily determined for t* — ► oo, which are valid 
close to the G-FP, within the mean-field-like PAD exten- 



two asymptotic and pre-asymptotic branches of a com- 
plete crossover effect. Therefore, any Fp [t*, D p (£*)] 
must contain at least 9 parameters. 

Eight of them can be readily unequivocally determined: 
e P from Eq. <2SJ, Z± from Eq. (2EJ, A from Eq. ll27l . 
A m / from Eq. 1(361) and two {X,, Yi} P pairs from the four 
Eqs. O, C2J, ©, and (SJJl. Using S$ >2 (A, A m /) as 
an entry data in Eq. (Jj|, leads to the supplementary 



sion t* > C 



mf 
PAD- 



In Eq. i?>'2\\ . A m j is the mean field 
exponent of the first order term of "classical" corrections, 
and Z Pm y is the associated amplitude defined by the 
following equation 



'P,mf 
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Y, 



p.i 



i=l X P,i 



(33) 



As in the case of the Ising-like PAD description, we 
impose four conditions on the ad hoc mean function 
to describe the mean-field PAD using the original pairs 
{X u Yi}p S of the bounded functions. These four condi- is composed of one crossover sensor S Pfl) character- 



determination of S P1 (A, A m f) from Eq. l!23ll . Sub- 
sequently, the parameters attached to the intermediate 
part of the crossover need to be determined. 

The phenomenological forms of Eqs. (0, © and ll23"l) . 
with K — 3, introduce 12 parameters [the case with K — 
2 (10 parameters) is of limited interest due to the error- 
bar propagation which cannot be accounted for via only 
two {Xi, Yi} p pairs]. Therefore, for any fitting procedure 
which uses S P2 (A, A m /) as an entry (free) parameter, 

only one {X t , Yi} p pair, remains truly free in the (K = 3) 
product terms of the MR crossover function. 

The triad {S-z,Xi,Yi} p of the calculated parameters 



tions read as follows 
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(34) 
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(35) 



(36) 



(37) 



Eq. H36H fixes the mean value of the mean-field-like 
"confluent" exponent, unequivocally. The three Eqs. 
ll34ll . lj35)l . and ll37jl . added to Eq. l(28ll . impose the mean 



values of two pairs among the 3 pairs {Xi,Yi} p 's. 



3.3. Number and nature of the parameters 



The above analysis demonstrates that the theoretical 
mean function F P \V%D p (£*)] must satisfy 9 constraints 
[Eqs. O, ED to 12HI, and $Q to JH3] to reproduce the 



istic of the t*-location of the classical-to-critical con- 
fluent crossover (see Figure EJ, and one "amplitude- 
exponent pair" {Xi,Yi} p , proper to the shape of the 
strict crossover part between the two asymptotic PAD 
behaviors (we will illustrate this latter point below using 
Figure H in § 3.4). 

Let us look for the possible existence of a unique value 
S2 = cte [or Si = cte by virtue of Eq. il23H ]. whatever 
the selected property and the considered (homogeneous 
or non-homogeneous) phase of the system. 

The S Pi (A,A m f) parameters are only related to the 
two universal confluent exponents. However, the numer- 
ical values of A appear also conditioned by the error- 
bar propagation of asymptotic uncertainties provided by 
the theoretical estimations of the leading universal expo- 
nents close to the non-Gaussian fixed point. In Figure 
we have reported by crosses at constant — y 1 ^- as a 
function of t* , all the [max, min, and mixing] conditions 
S P2 t/F — 1, whatever P and ± states (extending then 
the previous results reported in Figure for the homo- 
geneous susceptibility case). 

Since the dispersion of the mixing conditions 

\/ S $,max, 2 S $, min ,2 t * = l is significatively lowered, Fig- 
ure 3 supports the possibility of choosing a unique value 
for £2. 

Moreover, anticipating the S2 estimation given in the 
following subsection, the resulting mean value 1S2 = 
22.9007, common to all the properties and all the states, 



provides a condition S^v^o = 1 m close agreement with 
the previous mixed conditions, as shown by the corre- 
sponding vertical line in Figure The departure of 
the non-homogenous heat capacity case is probably due 
to the cumulative effect of the error bar propagation of 
the uncertainties on the exponent and on the critical 
background amplitude estimations, related here to the 
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Figure 3: Characteristic values of the crossover condition 
Sp 2 (A, A m /) yW — 1, for all the "min" (upper crosses at 
same Amax value) and "max" (lower crosses at same A m ,: n 
value) crossover functions estimated in Ref. (l4 | (lin-log A;i* 
diagram). The P (property) color indexations are: green 
X* (+), red t (+), blue C* (+), pink m* (— ), dark green 
X* (— ), and dark blue C* (— ). The corresponding values 
for the mixing crossover functions with E — \ [see Eq. 
I|29[l ]. are indicated by the intermediate crosses (at same 
A mean = m "^ """ value). The linear segments show 



the agreement with the associated geometrical mean values 
\/Sp max 2&p min 2^* — 1 of the min and max crossover con- 
ditions. The vertical line illustrates the "universal" crossover 



condition 52 ytjj = 1 [where D 



i+A r 



[S 2 (A„ 



from Eq. (jj], which is selected in the present 
study, whatever the property (P) or the phase domain (±). 



non-zero value of the regular background term (equal 
to 3) below T c . We can then note that the above well- 
controlled origin of this "universal" model-parameter S2, 
partly compensates for the arbitrariness of its numerical 
value. 



3.4. Unique form of the mean confluent function 



For each property, the three associated functions, 

(**) 



proposed in reference [Tj] can be considered. The three 
corresponding residual functions, rp (£*), expressed in 
%, are calculated from reference to our mean func- 
tion noted F P [t*,D(t*)], such as r P (t*) = 100 x 

1 I . The results reported hereafter are 
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1 P.mix 


t* D ± 
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Fp[t',D(t-)l 

obtained from the minimization method of the residuals 






P 


i/e 


i/x* 


C* 


1 


e+ 


-0.6303875 


-1.2395935 


0.1088375 


2 


z+ 


2.121008 


3.709601 


1.719788 


3 


A 


0.50189 






4 


z£+ 


-5.81623 


-8.56347 


8.06569 


5 


Si 


22.8573 






6 


s 2 


22.9007 






7 


Xp t i 


40.0606 


29.1778 


36.6874 


8 


Yp,i 


-0.098968 


-0.178403 


0.220033 


9 


Xp$ 


11.93211 


11.7625 


3.23787 


10 


Yp,i 


-0.153912 


-0.282241 


-0.000133095 


11 


Xp,3 


1.902735 


2.05948 


2.84102 


12 


Yp, 3 


-0.00789505 


-0.0185424 


-0.00222489 


13 


X c 






-3.79829 


14 


ep.mf 


-0.5 


-1 





15 


Z P,mf 


1 


1 


3.79004 


16 


A m/ 


-0.5 






17 


& P,mf 


-0.0195196 


-0.0391128 


0.00517327 



Table I: Numerical values of the parameters of the mean 
crossover functions Fp [t*,D(t*)] [see Eqs. 10 and 10, with 
K — 3 and S P2 — £2], corresponding to the dimensionless 
correlation length £*, susceptibility x*\ and specific heat C* , 
in the homogeneous phase (T > T c ). The lines 1 to 4 corre- 
spond to the characteristic parameters of the Ising-like PAD 
description [see Eqs. 12411 to 12811 ]. while the lines 14 to 17 
correspond to the characteristic parameters of the mean field- 
like PAD description [see Eqs. 13211 to 137j . The universal val- 
ues of the interrelated crossover sensors Si and S2 [see Eqs. 
12311 ], are given in the lines 5 and 6. The values of the three 
amplitude-exponent pairs {Xi,Yi} p given in lines 7 to 12, 
have been determined by a carefull adjustement to the theo- 
retical mixing functions of Eq. I|29|l (with E — |) proposed in 
Ref. [l4 |. using specific constraints of Equations 12511 to 128t 
and H34fl to H37fl . The value of the critical background Xc of 
the specific heat given in line 13 is the mean value of the min 
and max values of parameters X§ of Ref. [lj| . 



related only to the mixing function with E = | [see Eq. 

ESI- 

In a first step, starting with the susceptibility in 
the homogeneous phase as a basic property (since 
K = 3 already), we have validated the deriva- 
tion of several mean functions using several entry tri- 
ads made of \/ S^ max2 S^ min 2 and one among the 



three pairs {X, 
{^ ^ 



i.max-; * i,max L 5 i^ 



Yi. 



i,mtn 1 x t,mtn 






or 



-^i^max^i^mim 



c+Y it 



for each i value be- 



tween 1 to K. The two remaining pairs {Xj^i,Yj-a} s 

were then calculated from Eqs. (ESI, CSJl, CH) an d (53 - 
This step was then repeated for all the properties with 
two main results: 



1. at least one solution with K 
properties and all the states; 



3 exists for all the 
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P 


m* (MR7) 


m* (MR67) 


i/x* 


C* 


1 


e~ 


0.3257845 


0.3257845 


-1.2395935 


0.1088375 


2 


Zp 


0.937528 


0.937528 


17.762821 


3.203771 


3 


A 


0.50189 


0.50189 






4 


l}f 


3.42538 


7.70712 


-40.4666 


6.69984 


5 


Si 


22.8573 


22.8573 






6 


s 2 


22.9007 


22.9007 






7 


Xp,i 


219.597 


124.274 


470.671 


116.515 


8 


Yp,i 


-0.0284374 


0.0206744 


-0.0332665 


-0.105953 


9 


Xp,2 


38.3772 


20.9461 


43.6468 


72.9532 


10 


Yp- 2 


0.224464 


0.20453 


-0.585082 


0.247902 


11 


Xp,z 


6.92804 


6.92808 


5.2264 


12.6739 


12 


Vp,3 


0.152404 


0.123226 


0.139161 


0.0757262 


13 


X c 








-3.79349 


14 


e P,mf 


0.5 


0.5 


-1 





15 


^P.mf 


Vg 


V6 


2 


6.79349 


16 


A m .f 


-0.5 


-0.5 






17 


Z P,m.f 


0.0277176 


0.0277175 


0.0131516 


0.00846373 



Table II: Same as Table I for the dimensionless order pa- 
rameter m*, susceptibility x*i and specific heat C* , in the 
non-homogeneous phase (T < T c ). The distinction be- 
tween the MR7 and MR 67 results accounts for difference 
in the numerical values of the first confluent amplitude 1} M 
of Eq. 12SL which provides the (central) universal values 

■%■ (MR7) = 0.45 and ■%■ (MR67) = -% (MR6) = 0.9 

(see Appendix B for details). 



2. the residuals are minimum for the triad 

\f ^ 'max ,2^ min ,2> \f ^i, max ^M, mini 2 [ 



The following complementary observations were also 
made: 

1. The significative values of residuals are in a t* -range 
where S P2 ^/t* ~ 1. 

2. The t* -range where our mean functions compare to 
the max (or min) functions, corresponds to t* > 
-Qjr- (in this domain, the residuals are of the same 

order of magnitude than those obtained using the 
(E = i) mixing functions and remains < 1CP 3 ). 

3. Each Ising-like PAD description by equation 12411 
agrees with the complete crossover function within 
an error-bar lower than 0.01% for t* < C p s ^, 
where Cp^D ls then defined by the value 



S 



P:2 



rising 
-PAD 



10" 



(38) 



In a second step, for the inverse susceptibility and the 
inverse correlation length in the homogeneous state, we 
have minimized the residuals by successive small vari- 
ations around the above mixing values of each consti- 
tutive parameter of the entry triad. We have then se- 
lected the best mean functions for each property, using 



the mean value S2 = 22.9007 of the two optimized values 
S Xt2 = 22.9321 and % 2 = 22.8693, associated to \ and 
£, respectively. The final step was to minimize all the 
residuals with S% — 22.9007 fixed, whatever the property 
or the domain. As a main consequence, the universal con- 
fluent crossover condition S^v 7 ^* = 1 is such that t* = t, 
with 



0' 



/ 



1 



(s*y 



= 2.10" 



(39) 



Accordingly, Eq. IJ38JI has a universal form whatever the 
property, leading to well-controlled PAD extension t* < 



rlsinq ., i 
C PAD> wlth 



plsing 
'-PAD 



10 



-3 



w 



: w x tn = 2. 10 



-6 



(40) 



where w = 10~ 3 is then a convenient parameter to ac- 
count for the intermediate crossover range between the 
two (Ising like and mean field like) PADs (see Figure 
for the susceptibility case as a typical example ). 

The corresponding parameters of the theoretical mean 
functions are given in Tabled (homogeneous state, At* > 
0) and Table ITT1 (non-homogeneous state, At* < 0). The 
mean crossover functions for the susceptibility, the cor- 
relation length and specific heat in the homogeneous do- 
main are within the above theoretical level of precision in 
the complete £*-range. For the order parameter, suscep- 
tibility and specific heat in the non-homogeneous domain 
we get a sufficient accuracy (±0.1%) to have agreement 
with the theoretical crossover shape in the intermediate 
t*-range. The MR67 order parameter case in the non- 
homogeneous domain is the object of a specific analysis 
(see Appendix B) related to the application restricted to 
the one component fluid subclass. Correspondingly, for a 
more detailed analysis on the level of precision, we have 
also reported in this Appendix B the residuals for all the 
properties (see Figures 151 to llOl). 

To evaluate the influence of the constant universal 
value S*2 = cte — 22.9007, we have made a compari- 
son between all the t* ep -values where the following local 
mean value 



e P,eff 



e P 



"-mf 



(41) 



of the effective exponent ep, e jf occurs for any max, min, 
and mixing crossover function [the effective exponent is 

given by the equation ep, e // = l — ^— ^ iJ - [32JJ. 

In Figure^] using a vertical arbitrary unit (a.u.) scale 
between e Pmm i(a.u.) = to e Pmax> i(a.u.) = 1, then 

with e P i(a.u.) = e Pmix i(a.u.) = g at t* = t* , permits 
to generalize the (confluent) A case reported in Figure 
03 We can observe that the conditions of Eq. 1)4111 match 
perfectly at t* = t* ep for mean crossover functions and 
mixing crossover functions. 

Figure His the second result supporting our S% = cte 
suggestion to construct the "universal" confluent func- 
tion. The corresponding condition Ai = — 7^- for 
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CM 
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CO 



)0JC 



u 



10- 



10- 



10- 



e p 



position (semi-log scale) of the effective mean 
— — — (in arbitrary units), between critical 
Upper, lower, and me- 



Figure 4: 

value e p i = — ' ' 
and classical asymptotic exponents 
dian circles are for ep max, min, and mixing values of ref- 
erence |l3 |. while median crosses are for mean values of the 
present mean crossover functions. See the legend of Figure 
3 for the P color indexation. The vertical (black) line in- 
dicates the ^-position for the universal confluent condition 
S-2\/tfj — 1. Note the significative differences in t* ep -positions 
for max, min, and mean values of the effective leading expo- 
nents for each property, compared to the "universal" effective 
confluent exponent case reported in Figure 



the effective confluent exponent, where S2 y^ = 1 
[Eq. (J39l) ]. demonstrates that the t* ep -values are read- 
ily "C3" in nature (see Figure EJ, since, either the con- 

ditions ^(i) 2 < t* ep < 2(^) 2 , 

A e// (**J-Ai <I A 



or the conditions 



Asyvnp 



Ai 



with Aas 



syrap — 

AorA m f, are satisfied. We note that the relative in- 
differences between the values of t* (— ) w ^ ( ^- J , and 

t* (+) rs 2 I 4- j , for the inverse susceptibility in the 

non-homogeneous (— ), and homogeneous (+), phases, is 
correctly reproduced, while the t* -similarity, t* (— ) jw 

'oW w Hi") °^ the corres P on ding values for the 
specific heat, is also recovered (compare with the Figure 



CONCLUSION 



1 of Ref. [14]) 



By construction, the mean theoretical crossover func- 
tions determined in the present work do not account 
for any uncertainty on the parameters which, in the 
preasymptotic critical domain (PAD), characterize the 
critical behavior of the functions considered (correla- 
tion length, susceptibility, specific heat and coexistence 
curve). The recourse to such mean functions is justified 
by the fact that the amplitude and propagation of the 
asymptotic error-bars have a large extension which can 
affect significantly the classical-to-critical crossover espe- 



cially in the intermediate region where most experimental 
data are reported. Actually, it is known that in such a 
situation, when the comparison with experimental data is 
made following the common way of trying to determine 
the leading and first confluent amplitudes in the Ising- 
like PAD, the theoretical and experimental uncertainties 
add to each other. Consequently, it is usual to fix the ex- 
ponents to their mean theoretical values in order to get 
some definite information on the correction amplitudes. 
If we use the min and max theoretical crossover functions 
very close to the critical point, we obtain non-satisfying 
information on the intermediate crossover because the 
physical parameters $ and ip, introduced by the basic 
Eqs. ltl2)l and lTL3l) are not sufficiently well defined. The 
mean functions constructed with controlled constraints 
which correctly account for the universal features proper 
to the Ising-like PAD, provide useful criteria to define 
the extension of the Ising-like PAD. Fitting the experi- 
mental data with our mean crossover functions can then 
be a convenient method to validate any phenomenolog- 
ical approach by predicting the values of the two-scale 
factors which characterize the asymptotic singular be- 
havior within the Ising-like PAD and finally to identify 
eventual subclass of universality. In a forthcoming paper 
[]j|, the importance of the mean crossover functions will 
be illustrated by providing a complete characterization 
of the one component fluid subclass, within an extended 
asymptotic range well beyond the Ising like PAD, recover- 
ing then the practical intermediate range which is usually 
described by more complex formulations JS3, LlJ| ■ 



Appendix A: THE ADJUSTABLE PARAMETERS 



To simplify our analysis, we select the three- 
dimensional Ising-like < &^ =3 (n = 1) model, with go ~ 
[length] 

Let us consider the theoretical crossover function of 
Eq. 1J5J for the (dimensionless) correlation length £* (£*) 
given in Table []] at h* =0asa typical example, writing 
then 



and [see Eq. 1JT2J1] 



^* = .90 X £e 



t* = ti x At* 



(Al) 



(A2) 



Now fitting the experimental correlation length 
( el p(Ar*), in the homogeneous phase At* > 0, gives 
access to the two adjustable parameters 50 and $. Conse- 
quently, due to the asymptotic validity of Eq. l|A2l) close 
to the critical point, is known the asymptotic two-term 
Wegner expansion of the experimental singular behavior 
for the correlation length, which reads as follows 



£ exp = £+ (At*)" 



1 



(At* 



(A3) 
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Figure 5: Residuals (%) for the inverse susceptibility in the homogeneous phase. The right part compares relative residuals 
to the theoretical fitting precision of 0.01%. Color indexation: green (dark green) : MR7 (MR6) min crossover function; blue 
(dark blue): MR7 (MR6) max crossover function; red (dark red): MR7 (MR6) mixing crossover function with E — |; Residuals 
for MR6 functions are not visible in the magnified rigth scale, except in the C3 range. Four characteristics impositions are 
indicated by vertical lines which correspond to the t* < C-pad ('- e - the Ising like PAD extension) (orange line), the S^y^ = 1 
condition (black line), the t* = t* condition [Eq. I|41fl ] (blue line), and t* > £pad ('- e - the mean field like PAD extension) 
(pink line), respectively (see text for details). The amplitude and propagation of the asymptotic theoretical error-bars for the 
non-gausssian limit t* — » are clearly evident in the left part. To compare the amplitude of the theoretical error-bar to the 
one of the experimental uncertainty, the squared point and its associated theoretical error-bar indicate "equivalence" between 
the model and the experimental situation encountered at the "lowest " temperature distance T — T c — 1 mK above the critical 
temperature of xenon (see Ref. [lfil | for details). The rigth part corresponds to the residual magnification at the scale of the 
allowed [lj| theoretical precision level (< 0.01%). 
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Figure 6: Same as Figure 4 for the inverse correlation length in the homogeneous phase. 



The associated two-term expansion of the theoretical 
function reads as follows 



tween Eqs. IIA3J1 and IJA4L gives, 



r = (z|) (t*y v i-z)< + (t*r + 



(A4) 



Using Eqs. i|AHA2Jl . a term-to-term identification be- 



and 



-Z 1 / 



g^ 1 = $Z+V 



(A5) 



(A6) 



Within the preasymptotic domain, the two adjustable 
parameters of the model at h* = are unequivocally 
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Figure 7: Same as Figure 4 for the specific heat in the homogenous phase. Compared to results of Figures 4 and 5, the 
significative increases in amplitude and propagation of the error-bars is such as the residuals attain the admitted theoretical 
level in the C3 domain. 
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Figure 8: Residuals (%) for the order parameter in the non-homogenous phase. See Appendix B for details between the MR7 
(upper part) and MR 67 (lower part) results. 



defined by Eqs. l|A5Jl and IIAfiJl , Two important remarks 
can be made: 

i) the nonuniversal scale factor fl is uniquely defined 
by the confluent corrections to scaling; 



ii) the nonuniversal inverse length go is proportional to 
the inverse of the leading amplitude of the actual corre- 
lation length. 

>From now on, the complete characterization of the 
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Figure 9: Same as Figure 4 for the susceptibility in the non-homogenous phase. Compared to the general trend of all the other 
max and min crossover functions, a curious "exact crossing" between min and max functions occurs in the present case around 
t* m 310~ 13 , as revealed by the converging residuals for t* < 10~ 5 . That provides the maximun amplitude for theoretical 
error-bars observed outside the PAD extension. The associated error-bar propagation is then significative in a large C3 domain. 
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Figure 10: Same as Figure 4 for the specific heat in the non-homogenous phase. Compared to the homogeneous phase (Figure 
6), we suspect that the non-zero value of the backgroung term contributes to increase amplitude and error-bar propagation 
over the C3 domain, leading to its over-extended range close to the Gaussian limit t* — ■> oo. 



system would follow from the determination of ip which 
requires the consideration of one suplementary singular 
property, such as the susceptibility (as shown in the sec- 
ond paper 15]). 

However, the inverse length go, determined from cor- 
relation length measurements, is not a "natural" scale 
for the actual system because a measurement of the cor- 
relation length would be a prerequisite condition to set 
the length scale unit. One may easily show (23 that 
when the thermodynamic description of a magnetic sys- 
tem is normalized per particle, then all lengths are mea- 
sured in units of the thermodynamic microscopic length 



L mag 



(V 



mag J 



a/d 



(where v mag is the particle volume at 



criticality). In this respect, 



*"mag 



plays a role similar to 



that of the lattice spacing aj S i ng of the Ising uncompress- 



ible solid. The natural approach is thus to first choose a 
microscopic length unit such as aj S i ng , which implicitly 
introduces the number of particles per lattice cell. There- 
fore, the value of 50 is obtained via the length ai S i ng from 
the determination of an intermediate parameter u* Isi 
[like in Eq. lfl5ll ] such as: 



.90 X aising 



(A7) 



As a conclusion, fitting experimental results with the 
dimensionless theoretical functions requires reference to 
one length scale unit [ai S i ng , or a mag ] and the energy 
scale unit /cbT c to reduce to dimensionless quantities 
the thermodynamic and correlation functions. The fit- 
ting results of correlation length and susceptibility mea- 
surements for example [15J, enable the determination of 
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three dimensionless numbers ■&, ip, and go x aj S i ng (or 
.9o x a-mag)- In such a situation, d and ip act as the 
two-scale factors characteristic of the asymptotic univer- 
sal features of the Ising-like system, while the product 
go x a j S i ng (or go X a mag ) insures that the extensivity of 
the total macroscopic system is correctly accounted for 
in units of the actual critical correlation length. 

However, to compare between two systems of the same 
universality class remains not easy in the absence of 
explicit thermodynamic definition of the (coupling con- 
stant) inverse length g . This exercice is left to the second 
paper [15]. 



Appendix B: THE COEXISTENCE CURVE AND 
ITS CONFLUENT CORRECTION ERROR-BAR 



As explained in Ref. [Tj|, the MR7 crossover functions 
involve a forced account of the Guida and Zinn-Justin 
estimates [23 of the universal combinations between the 
leading critical amplitudes. This has induced an overesti- 
mation of the uncertainty on the correction terms in MR7 
compared to the previous MR6 calculations (see section 
IIB2 of Ref. (lj|). In general this does not have an im- 
portant impact on the resulting mean crossover functions, 
except for the order parameter for which the MR7 deter- 
mination finally appears to have (relatively to MR6) a 
poor quality (one may appreciate this difference by look- 
ing at bottom of Figures 2 (MR6) and 3 (MR7) of Ref. 
|22l |. for example). 



To circumvent this non-satisfactory error-bar situa- 
tion, we have applied the procedure described in Sec- 
tion 3 to construct a "modified" crossover function (la- 
belled MR67) of the order parameter. Actually, our 
MR67 function incorporates the same mean value of 
the leading amplitude as the MR7 one, Zm (MR67) = 
Z M (MR7) = 0.937528, and combines the MR6 central 

value of the universal ratio -r^r (MR6) = 0.9, with the 

MR7 amplitude value of 1%+(MR7) = 8.5635, so that 
1} M (MR67) = 7.70712. Using such a practice, we have 
reported the difficulty of accounting for error bars in the 
MR7 calculations of the crossover functions on only one 
property in the non homogeneous range. In the forthcom- 
ing paper dedicated to the study of the one component 
fluid subclass [ijj, a detailed analysis of this particular 
choice will be presented. 

The corresponding numerical values of the parameters 
of Eq. (§J , needed to calculate the MR67 mean crossover 
function for the order parameter in non-homogeneous 
state (t* < 0), are given in Table ITTI 

For a reader interested in checking the level of preci- 
sion for the present mean crossover functions in the three 
Ising-like, intermediate, and mean-field like, i*-range, we 



Fp,< 



[**] 



1 



have reported the residuals 100 x I — -p rpy 

for all the properties (see Figures 151 to I10JI . The right 
part of each figures gives the appropriate magnified scale 
of the residuals to compare with the estimated precison 
(±0.01%) for the fitting theoretical functions. 
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